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QUADRATURE FORMULAS FOR SIMPLE AND DOUBLE 
LAYER LOGARITHMIC POTENTIALS 

ELNUR H. KHALILOV AND MEHPARA N. BAKHSHALIYEVA 


Abstract. Quadrature formulas for logarithmic potentials and normal 
derivative of simple layer logarithmic potential are constructed. Error 
estimations for the constructed quadrature formulas are given. 


1. Introduction 

One of the methods for solving boundary value problems for Laplace equation 
is a method of boundary integral equations (see [2]). In many cases it is impos¬ 
sible to find exact solution of integral equations of boundary value problems for 
Laplace equation, therefore there arises a need to have justification of collocation 
method for these integral equations, and for this one should first construct quad¬ 
rature formulas for simple and double layer logarithmic potentials. Note that in 
[3], a cubature formula for simple layer acoustic potential was constructed, and 
in [4], a cubature formula for double layer acoustic potential was constructed. 
Then, based on these cubature formulas, the justification of collocation method 
for integral equations of boundary value problems for the Helmholtz equation in 
three-dimensional space was presented in [1, 5, 6, 7]. However, quadrature formu¬ 
las for simple and double layer logarithmic potentials have not been constructed 
yet. This work is dedicated to this problem. 


2. Dividing the curve into “regular” elementary parts 

Suppose a simple closed Lyapunov curve L <Z B? is given by the parametric 
equation x (t) = (xi (t) , X 2 (t)) , t G [a, b]. It is known (see [8]) that 


mi = V(^'i it)f + (4 {t)f > 0 


and 


Ml = maxy {x'^ (t))^ -|- {x '2 (t))^ < -boo. 

tS[a,fe] ' 


Let’s divide the interval [a, 6] into n > 2Mi {h — a) /d equal parts: 
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tk = a + , k = 0,n, where d is a radius of a standard circle (see [9]). 

As control points we consider x (rfc), where = a + A: = 1, n. Let 

Lk = {x{t) : tk-i <t<tk}, 

L'k = {x (t) : 4-1 <t <Tk} 

and 


Ll = {x (t) :: Tk <t <tk} ,k = 1, n. 
Then, using the arc length formula, we obtain 


mi {b — a) 
n 


< mesLk < 


Ml (6 - a) 
n 




and 


2n 

mi {b — a) 


< mesL'l. < 


2n 

Ml {b - a) 


2n - “ 2n ’ 

where mesL denotes the length of the curve L. Besides, taking into account the 
inequalities (see [8]) 

Mq mesL'f, < \x (r^) — x (4-i)| < mesL'j. 


and 

Mq mesL'l < \x {tk) — x (rfc)| < mesL’^, 

where Mq G (0,1) denotes a constant that does not depend on k and n, we obtain 

(1) V/c G {1, 2,n} : rfc(n) ~ Rk{n), where 


Tk (n) = min {lx {tu) - x (4-i)| \x {tk) - x (rfc)|} 

and 

Rk{n) = max{|x(rfc) -x(4_i)| , |x(4) -x(rfc)|}; 

(2) V/c G {l,2,...,n} : iifc(n) < (i/2; 

(3) V/c, j G {1,2, ...,n} : Vj (n) ~ rk (n). 

In (1), a{n) ~ b{n) Ci < < C 2 where Ci and C 2 are positive con¬ 

stants that do not depend on n. It is clear that r {n) ~ R{n) ~ where 
R {n) = max Rk {n) and 

fc=l, n 

r {n) = min (n). Besides, it is not difficult to prove the following lemma. 

fc=l, n 

Lemma 2.1. There exist the constants Cq > 0 and Cj > 0 independent of n 
such that yk,j G {1, 2,..., n} , j k and Vy G Lj the following inequality holds: 

CQ\y-x (rfc)l < |x (tj) - X (rfc)l < C'l \y - x {Tk)\ ■ 

In the sequel, we will call this kind of division a division of the curve L into 
“regular” elementary parts. 
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3. Quadrature formula for simple layer logarithmic potential 


Consider simple layer logarithmic potential 

P{x) = j <^ix, y) p{y)dLy, x = {xi, X 2 ) G L, (3.1) 

L 

where L <Z B? is a. simple closed Lyapunov curve with the index 0 < a < 1, p{y) 
is a continuous function on the curve L, <h(x,y) is a fundamental solution of the 
Laplace equation Art = 0, i.e. 

y) = 77 - 111 ]—— :,x,y G R^,x^ y, 

ZTT \x — y\ 

and A is a Laplace operator. 

Denote by C (L) a space of all continuous functions on L with the norm 
IIpII = max Ip (x) I, and introduce a modulus of continuity of the form 

xgL 

t VI V V ^ n 

0 ) = dsup-, 0 > U, 

T>S ^ 

for the function (p(x) G C (L), where cd{ip, r) = max \p{x) — <p(y)| . 

\x-y\<T 

x,y£L 

n 

Divide the curve L into “regular” elementary parts L = |J L^. 

fc=i 

Theorem 3.1. The expression 

h — ” /- 

Pn {x (Tfc)) = —^ J]] $ (x (Tfe) , X (Tj)) Y (x; {Tj)f + (x'2 {Tj)f p (x (t,)) 

is a quadrature formula for the integral (3.1) at the control points x(rfc), 
k = 1, n, and the following estimates are true: 

ui^\P{x{Tk))-Pn{x{Tk))\<M {u {p,l/n) + \\p\\^^) if 0 < a < 1, 

/c=l, n 

max |P(x(rfc)) - Pn (a:(rfc))| < M (w (p, l/n) + IIpII^ if a = 1. 

/c=l, n 

Hereinafter M denotes a positive constant which can he different in different 
inequalities. 


Proof. It is not difficult to see that 

P (x (Tfc)) - Pn {x (Tk)) = y $ (x (Tfc) , y) p (y) dLy+ 


+ 



a^(D))) p{y)dLy + 



p(x {Tj)))dLy + 
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n 

+ £ / ^ {X {Tk) , X (Tj)) (x; {t)f + (x'2 

7=1 . ^ 




-\/W ('0'))^ + K p{x{Tj))dt. 


Denote the terms on the right-hand side of the last equality by /i”(x(rfc)), 
/12 (x (t/j)), /i 3 (x (rfc)) and /14 (x (r^)), respectively. 

Using the calculation formula for curvilinear integral, we find 


\K{x{ru))\< 


{Tk), y)\dLy < 


R(n) 


In' - dr < M \\p\\^ R{n) \lnR{n)\. 


< M 


Let y G Lj and j 7 ^ k. Then, by Lemma 2.1, we have 


1 


1^ [x [Tk) ,y)-^{x {Tk ), X (tj))| = — 


1 

< — 
“ 27r 


In 1 


In ( 1 -h 
X {Tj) - y\ 


\x {Tk) - y\ 


\x{Tk) -x{Tj)\ - \x{Tk) -y\ 
\x{Tk) - y\ 

R{n) 


< 


< M- 


X {Tk) - y\' 


Consequently, 

\h^{x{Tk))\ < M \\p\\^R{n) J ^ < M \\p\\^R{n) \lnR{n)\ 
As the integral 


diamL 


r(n) 


|$(x, y)\dLy 


converges as an improper integral and 

J |< 1 > (x, 7 /)| dLy < M, Vx G L, 


(3.2) 


in view of Lemma 2.1 we obtain 

\h 3 {x {Tk))\ < MOJ {p, R {n)) J \^{x{Tk), y)\ dLy < Mu: {p,R{n)). 
Obviously, 

yj{x'l {t)f + {x'2 (t))^ - \l{x'^ (Tj))^ (x '2 {Tj))‘ 

< M {R{n))^ yt G [tj-i,tj]. 


< 


(3.3) 
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Then, taking into account the inequality (3.2) and Lemma 2.1, we find 


\K{x{t,))\< 



^>(x(rfc), x{Tj))\ dt< 


< 



(Tfc), y)\dLy < 


< M\\p\\^{R{n)f jj<l>{x (Tfc), y)\dLy < M ||/9||^ (i? (n))“ . 

As a result, summing up the estimates obtained above for the expressions 
hi{x{Tk)), /i .2 (ufc)), /13 (x(rfc)) and h2{x{Tk)), and taking into account the 

relation R (n) ~ we get the validity of Theorem 3.1. □ 


4. Quadrature formula for double layer logarithmic potential 


Now let’s construct the quadrature formula for double layer logarithmic po¬ 
tential 

W{x) = J p{y)dLy, x£L, (4.1) 

L 


where L C is a simple closed Lyapunov curve with the index 0 < a < 1, n{y) 
is an outer unit normal at the point y ^ L, p{y) is a continuous function on the 
curve L, and <k(x,y) is a fundamental solution of the Laplace equation. Divide 

n 


the curve L into “regular” elementary parts L 
The following theorem is true. 


U Lk. 

k=l 


Theorem 4.1. The expression 
b — a (a 

Wn (X (Tfc)) = — ^ -g 

i=i 

is a quadrature formula for the integral (4.1) at the control points x (r^), 


k = 1, n, and the following estimate is true: 

( In 71 \ 
Lu {p, 1/n) + IIpII^ — . 

71 J 

Proof. It is not difficult to see that 

^ (X {Xk)) - Wn (X (Tfc)) = j P ( 2 ^) + 

Lk 


(Xk) , x{Tj)) 
i{x (Tj)) 


>1 (d))^ + (4 {Xj)f P {X (Tj)) 



d4> (x (Tfc), x{Tj)) 
dn (x {Tj)) 


{p (y) - p{x {Tj)))dLy + 
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+E 




(x (Tk) , y) {x (rfc), x (tj)) 


dn{y) 


dn {x (tj)) 


P (y) dLy+ 


+E 


{x {Tk) , x{Tj)) 
dn {x {jj)) 




j-i 


X K (^)f + (4 (i)f - \l K + (4 p{x{Tj))dt. 

Denote the terms on the right-hand side of the last equality by (5”(x(rfc)), 
^2 (x (Tfc)) and 62 (x (rfc)), respectively. 

It is known that 


a4> (x (Tfc), y) 

1 

(yx{Tk),n{y)^ 

M 

5n (y) 


\x{Tk) -yf 

~ \x{Tk) - y|^““ 


Then, by the calculation formula for curvilinear integral, we obtain 

R(n) 


(x(Tk))\<M 


dr 


jz^<M\\p\\^{R{n)r. 


(4.2) 


Taking into account the inequality 


54> (x (Tfc) , y) 


dn (y) 


dLy < M 


dL't, 


\x{Tk) -y\ 


<M, V/c G {1,2, ...,n} 


and Lemma 2.1, we find 


\d 2 ix{rk))\ < Mu{p,R{n)) ^ 

j^k 

I '94> (x (Tk) , y) 


94> (x (Tfc), X (tj)) 
dn (x (tj )) 


dLy < 


< M uj {p, R (n)) 


dn{y) 


dLy < Moj {p, R (n)) 


Let y G Lj and j ^ k. It is clear that 


54> (x {Tk) , y) <94> (x {Tk) , X (rj)) 


9n (y) 


dn (x (tj)) 


+ 


+ - 


yx {Tk],n{y)^ (|a:(rfc) -x(tj)P - |x(rfc) - y|^) 
27r|x(rfc) -yp|x(rfc) -x(tj)|^ 

n{y)-n (x (r,))^ 
27r |X {Tk) - X {Tj)f 


+ 


27r |x (rfe) - X (rj)| 

Then from Lemma 2.1 we obtain 

I d<^ (x {Tk ), y) _ (94> (x (Tfc), X (rj)) 
dn{y) dn{x{Tj)) 


< 
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Consequently, 


R{n) {R{n)y^" {R{n)r 

\x{Tk)-yf~°' \x{Tk)-yf \x{Tk)-y\ 


1^3 ix{rk))\ < M IIpII^ X 


f atamL, 


diamL diamL \ 

/ %^mn)r f 


<AnlplL(H(n))“ |ln/;(n)|. 

Besides, from the inequalities (3.2), (3.3), (4.2) and Lemma 2.1 we have 


l«U(M)l<Af||p||^ (/;(»))“ 


1 


(9«h (X (Tfc) , X {Tj)) 
dn{x{Tj)) 


<A'/|ipiL(-R(")rE 


j^k 


d<^ (x (Tfc), y) 
dn (y) 


dLy < 


<M\\p\\^{R{n)r 


94» (x (Tfc), y) 
dn (y) 


dL,<M||y||^(i?(n))^ 


Summing up the estimates obtained above for {x (r^)), 5^ {x (r^)), 6^ (x (xk)), 
52 {x (ufc)) and taking into account the relation R (n) ~ i, we finish the proof of 
Theorem 4.1. □ 


Similarly, we can prove the following theorem: 

Theorem 4.2. The expression 

Qn (x (Tfc)) = ^dnlxiTk))^^^^ V W i^j)f + (4 P (Tj)) 

is a quadrature formula for the normal derivative of simple layer logarithmic 
potential 

„ . . dPix) f d^(x, y) . ^ 

^=A ~ahr ’ 

at the control points x{Tk), k = I, n, and the following estimate is true: 


max |Q(x(rfc)) - Qn (a;('rfc))| <M [uj{p, 1/n) + ||y||^ — . 

k=l,n \ 71 y 
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